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Abstract 

. The Global Colour Model(GCM) of QCD is used to extract low energy quark- 

CNl , gluon processes from experimental data. The resultant effective quark-quark 

coupling correlator is compared with that of Jain and Munczek, and with 
' the combined lattice results of Marenzoni et al. and Skullerud, and with the 

■ two-loop form. The results suggest that higher order gluon vertices are playing 

a role in coupling quark currents. The success of the GCM is explained by an 
infrared saturation effect. 
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1 Introduction 

We report on the current status of the project extracting low energy quark-gluon processes from 
experimental data using the Global Colour Model|J (GCM) of QCD. A new GCM fit (GCM98) 
to more extensive experimental data is reported here; this updates previous fits GCM95 and 
GCM97. The effective quark-quark coupling correlator is compared with that of Jain and 
Munczek ^j, with the combined lattice results of Marenzoni et a/.|| and Skullerud (J], and with 
the two-loop form. The results show that this correlator agrees remarkably well with that of 
Jain and Munczek, and with one constructed from the Marenzoni et al. and Skullerud lattice 
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results down to s — 1.8GeV 2 . But, significantly, all three of these depart from the two-loop 
form below s — 2.5GeV 2 . The difference between GCM98 - Jain-Munczek and the lattice 
construction below s = 1.8GcV 2 could be an indication of contributions to the quark-quark 
coupling at low energy from n > 4 gluon vertices. 



2 Global Colour Model 

The GCM modelling of QCD is based on the idea that as hadronic correlators are given by 
explicit functional integrals it should be possible, after an appropriate change of variables of 
integration, to identify a dominant configuration. It turns out that this dominant configuration 
is nothing more than the constituent quark(CQ) effect allowing QCD to be directly related to 
low energy hadronic physics. There are a number of reviews of the GCM|, |@, |,§. In the 
functional integral approach correlators are defined by 

= J DqDqDADCDC^.qjx) eM-SQCD{A^q,C,C]) 

jDqDqDADCDCexp(-S Q cD[A,q,q,C,C}) 

The various complete (denoted by scripted symbols) correlators Q lead to experimental observ- 
ables. They are related by an infinite set of coupled Dyson- Schwinger Equations (DSE), and 
by the Slavnov- Taylor gauge-symmetry related identities and, in the chiral limit, to the axial 
Ward-Takahashi identity (AWTI). The usual truncation of these DSE causes the violation of 
all these identities, and in particular of the AWTI leading to spurious effects for the dynamical 
breaking of chiral symmetry which is critical to low energy hadronic physics. One remedy for 
the AWTI problem is to modify the quark DSE so that, in conjunction with the meson equa- 
tions, the AWTI is satisfied. The GCM avoids this AWTI problem as it does not derive from 
these DSE; rather its nature follows from an analytical continuum estimation procedure for the 
functional integrations in which the AWTI is automatically satisfied. The correlators in Eq. (Q) 
may be extracted from the generating functional of QCD, 

Zqcd% V, J] = J DqDqDADCDCcxp(-S QC D[A,q, q,C, C]+rjq + QV + J A). (2) 

Functional transformations leadj^] to the GCM; briefly and not showing source terms for con- 
venience, the gluon and ghost integrations are formally performed 
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DqDqDADCDCcM~S QC D) = J DqDqcxp(- J q(-j.d + M)q+ 



+ f / J»#(I/)£W» -y) + §J H&fpGfvp + ). ( 3 ) 



where (x) = q(x) \-l^,q{x) , go is the bare coupling constant, and Q^ v (x) is the gluon correlator 
with no quark loops but including ghosts (C, C) 

_ = JPADCDCAl(x)At(y)exp(-S Q cD[A,C,C\) 
M J DADCDCcx P (-Sq C d[A,C,C}) ' [ ' 

Fig.l shows successive terms in Eq. ([}]). The terms of higher order than the term quartic in the 
quark fields are difficult to explicitly retain in any analysis. However the GCM models the effect 
of higher order terms by replacing the coupling constant go by a momentum dependent quark- 
gluon coupling g(s), and neglecting terms like Q^ p and higher order in Eq.(||). This g(s) is a 
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restricted form of vertex function. The modification g^G^ip) — > D^(p) — g(jp 2 )Q \iv(p)g(p 2 ) 
and the truncation then defines the GCM. We also call this effective quark-quark coupling 
correlator D M „(p) the effective gluon correlator. 






(a) (b) (c) 

Figure 1 : Diagrammatic representation of successive terms in the quark action in (J3|) . The quark-gluon 
vertex strength is go, while the gluon-gluon vertices (including gluon correlators) are fully dressed except 
for quark loops. 





(a) (b) (c) 

Figure 2: (a) The GCM effective _D M „ in (JE]), (b) example of correlations formally included in (a), and 
in (c) an n — 4 process not formally included in (a), but which is modelled in the GCM via the specific 
form of D M „. 

The GCM is equivalent to using a quark-gluon field theory with the action 

S GCM [A, q, q]= J (q(- r d + M+ iA^yJq + ^ApD~*(id)A* \ . (5) 

Here (p) is the matrix inverse of the Fourier transform of D^ v (x) and Fig. 2 shows processes 
formally included in D^^p). This action is invariant under q — > Uq,q — > qU\ and A^X a — > 
£/j4^A a [/^ (where U is a global 3x3 unitary colour matrix) - the global colour symmetry of 
the GCM. The gluon self-interactions that arise as a consequence of the local colour symmetry 
in Eq.(Q) and the ghost and vertex effects lead to D~^(p) being non-quadratic. Hence, in 
effect, the GCM models the QCD local gluonic action J F£ V [A]F£ V [A] in S QC d of Eq.(0) which 
has local colour symmetry, by a highly nonlocal action in the last term of Eq.(||) which has 
global colour symmetry. It is important to appreciate that while the GCM has a formal global 
colour symmetry the detailed dynamical consequences of the local colour symmetry of QCD 
are modelled by the detailed form of D(s). Approximations to the truncated DSE usually map 
QCD onto the GCM, as indicated in Fig.||. The success of this GCM modelling has been amply 
demonstrated §, |, @, |, 0. 

Hadronisation[|| of the functional integrations in Eq.(|I|) involves a sequence of changes of 
variables involving, in part, the transformation to bilocal boson fields, and then to the usual 
local hadron fields (sources not shown): 

Zw ( DqDqDAexp(SGCu[A,q,q\+fjq + qv) (GCM) 
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J DBDDDDcxp(-S[B, T>, £>]) (bilocal fields) 



= DNDN..DwDpDuj...exp(-S had [N,N,..ir,p,oj....}) (local fields) . 



(6) 



(7) 

The bilocal fields in Eq.(^J) arise naturally and correspond to the fact that, for instance, mesons 
are extended states. This hadronisation derives from functional integral calculus changes of 
variables which are induced by generalised Fierz transformations that emerge from the colour, 
spin and flavour structure of QCD. The final functional integrations in Eq.(|^) over the hadrons 
give the hadronic observables, and are equivalent to dressing each constituent hadron by, mainly, 
lighter constituent mesons. The basic insight is that the quark-gluon dynamics, in Eq. (jl|) , is 
fluctuation dominated, whereas the hadronic functional integrations in Eq.(Q) arc not. The 
induced hadronic effective action in Eq.(^|) is nonlocal. The GCM automatically preserves the 
consequences of the dynamically broken chiral symmetry in the action in Eq.(^). 



Truncated SDE 




QCD 



GCM 



NJL, ChPT 



Hadronisation *- OBSERVABLES 




MIT, Cloudy Bag 
Solitons, QHD, QMC 



Lattice »- Hadron Correlations' 

Figure 3: Relational map of the GCM to QCD and various other modellings including the Nambu 
- Jona-Lasinio (NJL), Quantum Hadrodynamics (QHD), Quark Meson Coupling (QMC) and Chiral 
Perturbation Theory (ChPT). 

A key idea in the GCM is that in proceeding from Eq.(^|) to Eq.(0) one expands S[B, ..] 
about the configuration Bcq that minimises it, giving in Eq.(||) the GCM Constituent Quark 
(CQ) equations, 
SS , 



Thus for all hadrons one assumes a universal dominant configuration. This amounts to assuming 
that all hadrons share a common dynamical feature. Of the set B(x,y)cQ only A(x — y) 
and B(x ~ y) (their Fourier transforms appear in Eq.([lO|)) are non-zero translation- invariant 
bilocal fields characterising the dominant configuration. This is the dynamical breaking of chiral 
symmetry. The dominant configuration is analogous to the condensate of Cooper pairs in the 
BCS theory of superconductivity. However unlike solid state superconductivity here there is 
no normal component. Writing out the translation invariant CQ equations we find that the 
dominant configuration is indeed simply the constituent quark effect as they may be written in 
the form, 

4 r d 4 q 
3 J J2nj' 

and we see that this is the gluon dressing of a constituent quark. This equation is exact in the 
GCM and its solution has the structure 



G 1 {p) =ifr + m + - 



(9) 



G{q) = {iA{q)q. 1 + B(q) + m) 1 = -iq.-f<j v (q) + a s {q). 



(10) 
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In the chiral limit there are more Bqq fields that are non-zero, and the resultant degeneracy 
of the dominant configuration is responsible for the masslessness of the pion. The constituent 
quark correlator G should not be confused with the complete quark correlator Q from Eq. ([!]). 
This Q would be needed to analyse the existence or otherwise of free quarks. The G on the 
other hand relates exclusively to the internal structure of hadrons, and to the fact that this 
appears to be dominated by the constituent quark effect. The evaluation of Q is a very difficult 
task, even within the GCM while G is reasonably easy to study using Eq.(||). The truncation of 
the DSE in which the full quark Q is approximated by the CQ G amounts to using a mean field 
approximation; however the truncated DSE then has no systematic formalism for going beyond 
the mean field as in the GCM. The hadronic effective action in Eq.(^J) arises when S[B, ..] is 
expanded about the dominant CQ configuration; the first derivative is zero by Eq.(|^), and the 
second derivatives, or curvatures, give the constituent or core meson correlators G m (q,p; P) 

G m 1 (q,p;P)=F.T. (—-^- \ B ) , (11) 



6B(x,y)SB(u,v) ]BcQ 

after exploiting translation invariance and Fourier transforming. Higher order derivatives lead 
to couplings between the meson cores. The G m (q,p; P) are given by ladder-type correlator 
equations. The non-ladder effects are inserted by the final functional integrals in Eq.(^), giving 
the complete GCM meson correlators Q m (q,p; P). It is interesting to note that the truncated 
and modified DSE in Maris and Roberts p3[ are identical to Eq.(pl) and Eq.(|ll|) (in the form 



of Eq.(12)). However Maris and Roberts |13|] assume that D(s) has the two loop form down to 
s w lGeV 2 , and this assumption is not supported by the detailed analyses reported herein. An 
important advantage of the GCM modelling of QCD is that the hadronic sector effective action 
is manifestly derivable. One key feature is that the meson exchange within the diquark subcor- 
relations of the baryon correlators is seen to insert crossed gluonic exchanges, and these appear 
to stop the diquark correlations from developing a mass-shell, i.e. the diquark correlations are 
confined fl4||. 

In the present analysis the u and ai mesons are described by these constituent meson 
correlators; that is, we ignore meson dressings of these mesons. The mass M of these states 
is determined by finding the pole position of G m (q,p; P) in the meson momentum P 2 = —M 
and this leads to the homogeneous vertex equation 

r(p; P) = ~\f - v)iM<i + f )r(g; P)G( q - |) 7 ,. (12) 



The success of the GCM appears to be based on the phenomenon of an IR saturation 
mechanism in which the extreme IR strength of the many contributing quark-quark couplings 
is easily modelled by this one effective gluon correlator D^ip). Of particular dynamical interest 
is the comparison of the GCM D fiU {p) with one constructed theoretically from only a gluon 
correlator and vertex functions, say from continuum or lattice modellings, for this gives some 
insight into the IR strength of the higher order gluonic couplings. 

3 Procedures 

To solve Eq.(9h for various D llv {p) and then to proceed to use A[s) and B(s) in meson correlator 
equations for fitting observables to meson data is particularly difficult. A robust numerical 
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technique is to use a separable expansion [[l0|. 

D^(p) = (S^-^)D(p 2 ), and ^(p) = (6^ - ^)G(p 2 ). (13) 

We expand D(p — q) in Eq.(^) into 0(4) hyperspherical harmonics 

D{p -q) = D (p 2 , q 2 ) + q.pD.ip 2 , q 2 ) + (14) 

D {p 2 ,q 2 ) = - [ df3sm 2 (3D(p 2 + q 2 - 2pqcQS0), ... (15) 
n Jo 

We then introduce multi-rank separable Dq expansion (here n = 3): 

D ( P 2 ,q 2 )= J2 ?i(p 2 )ri(q 2 )- (16) 

i=l,n 

The constituent quark equations then have solutions of the form 

B(s) =J2 B i(s), B((«) = 6 4 ri(a), a a (s) = ]T ^(s);, (17) 

i—l,n 

b 2 = Air 2 I™ sdsB^asis), and B^a) = (18) 

However rather than specifying Ti in (|l^) we proceed by parametrising forms for the cr SJ ; 
and cr v ; the then follow from ( p^ ) and ([l8|): 

^ 2.s 2 -rf 2 (l-exp(-2.s 2 /rf 2 )) 
ct s (s)i = ciexp(-dis), (J s (s) 2 = c 2 . I ^"-j 

, ^ f 2/(s)-d 3 (l-exp(-2/( 8 )/d 3 )) V , , , /A2 ^ 1/2 

ct s (s) 3 = c 3 I 2f(sf ) ' ' ' ^ = " ' ' ' 

2 5 -/3 2 (l-exp(-2 s //3 2 )) 
cr„(s) = — 2 . (19) 

The three a S i terms mainly determine the IR, midrange and UV regions; the cr s (s) 3 term 
describes the asymptotic form of a s (s) ~ l/s 2 ln(s/A 2 ) for s — * oo and ensures the form for 
-D(s) ~ 1/s ln(s/A 2 ). With these parametrised forms we can numerically relate the mass of the 
ai(1230MeV) and w(783MeV) mesons from Eq.@, U(93.3MeV) and experimental points for 
a(s) (see Fig. 4a insert) from the Particle Data Book for s > 3GeV 2 to the parameter set in 
Table 1 in a robust and stable manner. 

Table 1: a s (s) and cr v (s) Parameters. 





1.839GCV- 1 


di 


3.620GeW 2 


P 


0.4956GeV 


c 2 


0.0281GeV 7 




1.516GeV 4 


A 


0.234GeV 


C3 


0.0565GeV 3 


d 3 


0.791 lGeV 2 







The translation invariant form for the effective gluon correlator is easily reconstructed by 
using D(p 2 ) = D {p 2 ,0) from Eq.@ 

V b 2 ct s (0) 2 p 2 a v (p 2 ) 2 + <J s {p 2 ) 2 



With the parameter set in Table 1 the resulting D(p 2 ) is shown in Fig.[|. Shown in Fig.^a is the 
pure gluon correlator Q{p 2 ) from the lattice calculations of Marenzoni et al [pi, corresponding 
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to the value /3 = 6.0, where the errors arise mainly from a 5% uncertainty in the lattice spacing; 
a = 0.50±0.025GeV _1 . A significant feature of QCD is that the infrared dominance, as revealed 
by the large value of D(s) at small s, causes the CQ equations to saturate, i.e. the form of the 
solutions A(s) and B(s) is independent of the detailed IR form of D(s). This saturation effect 
means that low energy QCD is surprisingly easy to model. 





(a) s(GeV 2 ) (b) s(GeV 2 ) 

Figure 4: Plots of D(s) (GeV~ 2 ). In (a) solid line is GCM98; shortdash line is GCM95; longdash 
line is GCM97. Data plot is lattice pure-gluon Q(s) from Marenzoni et al., and so has no quark- 
gluon vertex. Insert is fit of GCM98 to a(s) of Particle Data Book. In (b) GCM98 is the solid 
line; dashed line is Jain and Munczek; lower solid line is two-loop form with A = 0.234GeV and 
Nf — 3; data plot is combined lattice data for g(s)G{s)g(s) with G(s) from Marenzoni et al. 
(as in (a)) and g(s) from Skullerud. Insert shows <?(s) from Skullerud (lower data plot), and 
from GCM98/£(s) (upper data plot). 

Also shown in Fig.|(a) are the earlier GCM fits: GCM950 and GCM970. The insert in 
Fig.|](a) shows the GCM98 in the form of a(q) {D(s) = inHifi) compared with experimental 
data from the Particle Data Book. In Fig.[|(b) we show the new GCM98 quark-quark coupling 
correlator D(s) showing excellent agreement with the Jain-MunczekQ D(s), and with one 
constructed from the Marenzoni et al. and Skullerud lattice results down to s = 1.8GeV 2 ; 
D(s) = g(s)G(s)g(s). The normalisation and shape of the Marenzoni et al. G(s) agrees with 
that of Leinweber et a/.|15|. However the normalisation of the much more difficult lattice 



computation of g(s) is uncertain and we have chosen it so that the combined lattice D(s) 
agrees with the experimental Particle Data Book for s > 3GeV 2 . As shown in Fig^(b) all 
three of these depart from the two-loop form below s = 2.5GeV 2 . The difference between 
GCM98 - Jain-Munczek and the lattice construction could be an indication of contributions to 
the quark-quark coupling from n > 4 gluon self-couplings at low energy, since processes like 
Fig.||(c) would be included in the GCM fit, but are not in the lattice construction. 
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The insert in Fig.^(b) shows the g(s) from g 2 (s) — D(s)/Q{s) that then follows from our 
analysis. This is the effective quark-gluon coupling vertex if the gluon correlator is taken to 
be that of Marenzoni et al. or Leinweber et at. Here the error bars now indicate combined 
errors and uncertainties from the lattice spacing. Also shown is g(s) from SkullerudQ with the 
normalisation as discussed above. 

The GCM95-GCM98 parametrisations differed mainly in their asymptotic forms, and when 
fitting to meson data these differences only resulted in slight variations of D(s) for intermediate 
s values, as shown in Fig.|](a). The consequent small variations in the predicted values of 
various hadronic properties are shown in Table 2. However forcing the asymptotic form to fit 
experimental data for D(s) even for s > 3GeV 2 results in a stabilisation of the GCM D(s). 



Table 2: Hadronic Observables. 



Observable 


GCM1995 


GCM1997 


GCM1998 


Expt. /Theory 


A 


93.0MeV* 


93.2MeV* 


92.40MeV* 


93.3MeV 


a\ meson mass 


1230MeV* 


1231MeV* 


1239MeV* 


1230MeV 


7r meson mass (for m U] d) 


138.5MeV* 


138.5MeV* 


138.5MeV* 


138.5MeV 


a(q) 


~ 


- 


see fig. 4a 


t 


K meson mass (for m s only) 


496MeV* 


- 


- 


496MeV 


(m u + m d )/2\ R (Li= lGeV) 


6.5MeV 


4.8MeV 


7.7MeV 


R;8.0MeV 


m 8 \ R (n = lGeV) 


135MeV 






130MeV 


lu meson mass 


804MeV 


783McV* 


783MeV* 


782MeV 


Oq it — it scatt. length 


0.1634 


0.1622 


0.1657 


0.26 ± 0.05 


do it — 7r scatt. length 


-0.0466 


-0.0463 


-0.0465 


-0.028 ± 0.012 


a\ n — 7r scatt. length 


0.0358 


0.0355 


0.0357 


0.038 ± 0.002 


a-2 it — it scatt. length 


0.0017 


0.0016 


0.0018 


0.0017 ± 0.003 


a! it — it scatt. length 


-0.0005 


-0.0005 


-0.0003 


0.00013±0.0003 


r-n pion charge radius 


0.55fm 


0.53fm 


0.53fm 


0.66fm 


| + (0 + )nucleon-core mass** 


1390MeV 


1435MeV 


1450MeV 


>1300MeVtt 


constituent quark rms size 


0.51fm 


0.39fm 


0.58fm 




chiral quark constituent mass 


270MeV 


267MeV 


325MeV 




+ diquark rms size 


0.55fm 


0.55fm 


0.59fm 




+ diquark const, mass 


692MeV 


698MeV 


673MeV 


>400MeV 


1 + diquark const, mass 


1022MeV 


903MeV 


933MeV 




CP diquark const, mass 


1079MeV 


1049MeV 


1072MeV 




1~ diquark const, mass 


1369MeV 


1340MeV 


1373MeV 




MIT bag constant 


(154MeV) 4 


(145MeV) 4 


(166MeV) 4 


(146MeV) 4 


MIT N-core (no cms corr.) 


1500MeV 


1420MeV 


1625MeV 


>1300MeVft 


* fitted observable; - not computed or not known; t q(s) from Particle Data Book; 


GCM1995: Cahill and Gunnerjli 


1; GCM1997: Cahill and GunnerO; GCM1998: this report. 


** only + diquark correlation; 1 


+ diquark correlation lowers nucleon core 


mass. 


ft nucleon core mass (i.e. no meson dressing) 









4 Conclusion 

Figs|4](a) and |^(b) reveal that the GCM95-98 project of extracting low energy quark-gluon 
processes has reached a stage where some detailed insights are emerging. It is clear that the 
continuum modelling, when fitted to experimental data, leads to a unique quark-quark coupling 
correlator. This is shown by the excellent agreement between the GCM multi-rank technique 
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and the Jain-Munczek result [£| O. This continuum result is also in excellent agreement with 
the combined lattice prediction down to s rj 1.8GeV 2 . All three results show a clear departure 
from the two-loop form below s — 2.5GeV 2 . The difference between the GCM modelling and 
the combined lattice D(s) at lower s values raises an interesting interpretation. Assuming that 
the lattice calculation of g(s) is confirmed by future studies, particularly at low s, we must 
conclude that the GCM analysis of the experimental data is revealing the influence of processes 
that are not in the lattice construct. The most obvious possibility is that the high order gluon 
vertices, as in Fig.||(c), are contributing at low energies, and providing additional attractive 
interaction between quark currents. Provided these contribute only in the region where D(s) 
is large the detailed form of these contributions is not needed since the saturation effect means 
that hadronic observables are independent of these details, including their gauge dependence. 
Thus the saturation of the constituent quark processes and the related dynamical breaking of 
chiral symmetry explains why the GCM works so well as a low-energy equivalent field theory 
for QCD. 
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